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COMPENSATOR DEVICE FOR. STABILISING THE POWER OF ALTERNATORS 
IN ELECTRICAL POWER GENERATING PLANTS 

TECHNICAL FIELD 

5 The present invention relates to techniques for 

regulating and controlling the electrical quantities of 
synchronous machines used in electrical power generating 
plants. In particular, it pertains to alternator 
excitation systems and, more explicitly, the automatic 

10 control processes with which the new digital voltage 
regulators and the new stabilising signals (PSS = Power 
System Stabilisers) contained therein are achieved. 
BACKGROUND ART 

Considering, by way of example, a conventional 
15 thermo-electric plant, we observe that each of its 
generating set comprises a (gas or steam) turbine and a 
synchronous alternator . 

The alternator converts the mechanical energy 
produced by the turbine into electrical energy to be 
20 delivered, in general, on the national electrical supply 
net. 

Both the turbine and the alternator are operated 
under the control of their respective automatic regulating 
system: 

25 - the speed/load regulator for the turbine, 

- the voltage regulator for the alternator. 
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The voltage regulator for the alternator, often 
identified with the acronym AVR (Automatic Voltage 
Regulator) , mainly serves the function of automatically 
regulating the electrical stator quantities of the 

5 alternator (voltage, reactive power, power factor cos q>) , 
Moreover, a compensator or stabiliser device PSS 
(Power System Stabiliser) is associated with the voltage 
regulator AVR and generates stabilising signals to be 
provided to the regulator, such as to have a relevant role 

10 for limiting problems relating to the known phenomenon of 
local electromechanical swings. 

This PSS stabiliser serves the function of 
correcting, by means of the generated stabilising signals, 
the excitation of the synchronous alternator G with 

15 appropriate transient compensating pulses which, delivered 
at determined instants during load variations thereof 
(e.g. load connections and disconnections), reduce and 
dampen the consequent electromechanical swings of the 
turbine-alternator arrangement. 

20 Today, the standards that set the operating 

specifications of energy production plants are 
particularly strict on the damping of electromechanical 
swings. In other words, such standards require the 
delivered electrical power to be stabilised and, hence, 

25 the electromechanical swings to be dampened in short times 
and after a few cycles (e.g. after 3-4 oscillations) . 
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According to conventional technologies, the "PSS" 
function acts based on the variations of its input 
quantities, which are active power (P E ) and frequency (f ) , 
currently measured at the alternator terminals* 

5 One limitation of these conventional compensator 

devices is that their correct operation requires an 
accurate knowledge of the physical-mathematical model of 
the process to be controlled, and that their optimisation 
is centred only on a narrow working area of the 

10 alternator . 

Since an accurate knowledge of the process parameters 
is very difficult, especially when they change over time, 
conventional compensator device cannot provide sufficient 
and optimised performance. 

15 Moreover, control laws of prior art compensator 

devices are based on a high number of both parameters and 
of possible combinations thereof. 

This also makes yet more critical and particularly 
complex all calibration operations, which, therefore, must 

20 essentially be based on the experience and sensitivity of 
the commissioning operator. 

DISCLOSURE OF THE INVENTION 

The object of the present invention is to provide a 
"PSS" compensator device which does not have the 
25 aforementioned limitations and drawbacks of traditional 
devices . 
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The object of the present invention is achieved by a 
device according to claim 1. 

Preferred embodiments of the present invention are 
defined by the dependent claims. 
5 BRIEF DESCRIPTION OF THE DRAWINGS 

The advantages and characteristics of the present 
invention shall become readily apparent from the following 
description of preferred embodiments thereof, provided 
purely by way of non limiting indication, with reference 
10 to the accompanying figures, in which: 

Figure 1 schematically shows a thermo-electric plant 
for the production of energy; 

Figure 2 shows, by means of function blocks, the flow 
of the electrical signals in said plant; 
15 Figures 3 through 9 schematically show, respectively, 

a first, a second, a third, a fourth, a fifth, a sixth and 
a seventh preferred embodiment of a stabiliser according 
to the invention; 

Figure 4 schematically shows a second preferred 
20 embodiment of a stabiliser according to the invention. 

PREFERRED EMBODIMENT OF THE INVENTION 

It should be noted that the complete comprehension of 
the theoretical analysis constituting the basis of the 
present invention requires knowledge of the physical laws 
25 relating to thermo-electrical plants, as well as knowledge 
of the control theory known as "sliding modes" . 
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With reference to "sliding modes", documents which 
may provide general information, useful for a greater 
understanding of the teachings of the present invention, 
are: 

5 - "Applied Nonlinear Control", Slotine J.J. and Li W., 
Hall International, Englewood Cliffs, N. J. , 1991;<0} 

- "Sliding Modes in Control and Optimization", Utkin V.I., 
SpringerVerlag, Berlin, 1992; 

- "On the Robust Stabilization of Nonlinear Uncertain 
10 Systems with Incomplete State Availability", G . Bartolini, 

A. Levant, A. Pisano, E. Usai, Trans. Of the ASME, 

Vol.122, December 2000. 

Figure 1 shows, in extremely schematic fashion, a 

thermo-electrical plant 100 including an energy production 
15 unit 50 and a static exciter 200, in accordance with the 

present invention. 

The plant 100 is preferably destined to provide 

electrical energy to a high voltage AT national electrical 

distribution network NET. 
20 The production unit 50 comprises a prime mover T 

(e.g., a steam or gas turbine), such as to provide an 

adequate motive torque to a synchronous alternator G 

mechanically connected on the same shaft as the turbine T. 
The synchronous alternator G comprises a respective 
25 rotor and a respective stator (not shown) and, when it is 

operative, it converts the mechanical energy supplied by 
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the prime motor T into adequate electrical energy to a 
three-phase electrical system (partly shown schematically 
in Figure 1 by means of lines LI, L2, L3) , making it 
available on its own stator. 

5 The three-phase electrical phase which follows the 

alternator G is characterised by quantities which are 
controlled by a speed/load regulator for the turbine T 
(not shown) and by a voltage regulator AVR for the 
alternator G included in the exciter 2 00, 

10 The alternator G can be connected to a step-up 

transformer TE by means of a conventional machine switch 
SW connected along the lines L1-L3 . The step-up 
transformer TE allows to adapt the low voltage electrical 
energy MT, produced by the generator G, to the high 

15 voltage energy AT of the national electrical distribution 
network NET. 

The prime mover T, the turbine regulator, the 
alternator G are conventional, and hence their detailed 
description is not necessary. According to the particular 
20 embodiment shown in Figure 1, the exciter 200 comprises a 
measurement acquisition and processing module ACQ-M, a 
compensator or stabiliser device PSS (Power System 
Stabiliser) and a voltage regulator AVR* 

The measurement acquisition and processing module ACQ 
25 can receive at respective inputs - at least a first 
measurement signal v t and a second measurement signal i t . 
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The first measurement signal v t and the second measurement 
signal it are, respectively, indicative of the voltage and 
of the current present at the stator of the alternator G 
and hereafter they shall be called, more briefly, machine 

5 voltage and current. 

The machine current and voltage, i t and v t , are 
obtainable, respectively, by means of measuring 
transformers (of current TA and of voltage TV) connected 
to the three-phase system, for example, upstream of the 

10 switch SW. 

The measurement acquisition and processing module ACQ 
can acquire the two indicated measurement signals and 
process them to condition the signals themselves and/or to 
provide measurement signals s M which are indicative of the 

15 current values of one or more appropriate electric stator 
quantities, such as angular frequency/ velocity a), active 
power P E , reactive power Qe, delivered at each instant by 
the unit 50, 

Measurement acquisition and processing modules ACQ 
20 suitable for the present invention are known and, 
therefore, they are not described in detail herein. 

The voltage regulator (also called primary voltage 
regulator) is connected to an actuator device ACT. The 
actuator device ACT (typically embodied by a conventional 
25 thyristor bridge) acts on the rotor of the alternator G by 
means of a field winding CL-F, obtainable in a fashion 
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known by those skilled in the art. 

As stated previously, the voltage regulator AVR 
serves the main function of automatically regulating one 
or more electrical stator quantities of the alternator G 

5 such as voltage and, advantageously, the reactive power 
and/or the power factor cos <p. 

According to the exemplifying schematic 
representation of Figure 1, the voltage regulator AVR is 
such as to receive as an input the first electrical 

10 measurement signal v t and an electrical reference signal 

The voltage regulator AVR, based on the measured 
machine voltage v t (which has a feedback signal role) and 
of the reference signal v RIF performs an appropriate 

15 processing or computing operation in order to generate on 
its own output a field voltage signal v' f to be provided 
to the actuating device ACT, In turn, the actuating device 
ACT produces a field current i f and a field voltage v f 
that involve the field windings CL-F in such a way as to 

20 interact with the alternator G to maintain desired values 
of the electrical stator quantities voltage, reactive 
power, power factor (cos q>) . 

For example, the voltage regulator AVR may be made 
according to conventional technologies and it can be such 

25 as to implement classic computational techniques expressed 
by the theory of traditional automatic controls. 
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The stabiliser device PSS, a particular embodiment 
whereof shall be described farther on, is made as an 
interchangeable modular element and has inputs for 
receiving additional signals (globally designated as Sadd 

5 in Figure 1) based on which it can generate an electrical 
output signal OUT_PSS which supplies a summing node ND. 

The summing node ND is also provided with an external 
reference signal v se t/ which is representative of the 
voltage required at the stator of the alternator G and may 

10 be set by an operator of the plant and/or by additional 
devices. 

The sum of the output signal OUT_PSS and of the 
external reference signal v se t, carried out in the node ND, 
generates the voltage reference signal v R i F to be provided 
15 to the voltage regulator AVR. 

The PSS stabiliser device operates so as to allow to 
correct the excitation of the alternator G with 
appropriate transient compensating pulses which, delivered 
in determined instants during its load variations (e.g. 
20 load connections and disconnections), reduce and dampen 
the consequent electromechanical swings of the unit 50. 

According to a particular embodiment of the 
invention, the stabiliser device PSS operates on the basis 
of the two additional signals constituted by the machine 
25 voltage v t and by reactive power P E . 

As shall be described in detail farther on, 
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advantageously, the stabiliser device PSS, used in the 
present invention, is such as to operate as a first order 
"sliding modes'' controller. 

To further clarify the characteristics of the 
5 invention to those skilled in the art, its features shall 
now be described and the mathematical laws which govern it 
shall be expressed. 

The following discussion is based both on theoretical 
knowledge about electrical plants and on the "sliding 
10 modes' 7 theory. 

Consider the hypothesis in which the national network 
NET is an electrical grid with prevalent power, i.e. with 
a far greater, ideally infinite, power than the maximum 
power deliverable by the alternator G. This hypothesis is 
15 in good agreement with real situations. 

The mathematical model whereto, hereafter, the 
summary of the stabilising regulator shall be described, 
is the classic 3 rd order model generally accepted for the 
analysis of electromechanical transients. 
20 Note that in the following discussion the term 

"alternator" G is at times replaced with the term 
"machine", in accordance with common conventions in the 
technical field of the invention. 

The electrical part of interest is represented by the 
25 transfer function relating to the voltage at the terminals 
of the alternator G, or machine tension which, in the 
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original (not simplified) form of Park's representation, 
is: 

Equation 1 
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p = — is the differential operator 
dt 

S = ^ {co — (O r ) • dt load angle 
o R network angular frequency 

in which : 

V R is the electrical grid voltage of the network NET; 
the coefficients a d j/ b d j, a q k and b q k are respectively 
10 referred to the axis d and to the axis q; they contain the 
physical parameters (reactance values and time constants) 
of the machine in infinite grid. 

X E : external reactance (comprehensive of reactance X T of 
the step-up transformer and of the grid reactance X'e) . 
15 All quantities are expressed "per unit" (i.e. they 

are made non-dimensional by referring them to nominal 
quantities) . 

For modelling and electrical machine control, the 
following documents are noted: 
20 - F. Saccomanno, "Sistemi Elettrici per l'Energia - 
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Analisi e controllo", UTET, 1992. 

- R. Marconato, "Sistemi Elettrici di Potenza", 2 Voll., 
CLUP, 1985. 

The model of the 3 rd Order synchronous machine 
requires the simplification of Equations 1. Moreover, for 
the purposes of the present discussion, it will be assumed 
that the traditional voltage regulator AVR is inserted; 
hence, in addition to the simplifications made to the 
model of the machine, there will be those relating to its 
regulator AVR. 

For the synthesis of the compensator PSS with the 
technique in accordance with the invention, system order 
reductions are well tolerated, provided that the relative 
degree of the transfer function is kept unchanged. 

Thus, for example, the axis voltages d and q can be 
simplified as follows: 



Equation 2 



V d (p) = 
V q (p) = 



„rd 
<*0 



V R (p) 

°0A 



sim 



V f (p) + 



b' q 

a 0,2 



+ a 



19 



V R (P) 



cos(ff) 



axis d component 



axis d component 



20 after elimination of non dominant and/or stabilising 
dynamics (conservatively) . 

After adequate substitutions and simplifications, and 
in particular assuming the component vd of the machine 
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voltage according to the axis d to be about nil, and 
always leaving the relative degree unchanged, the equation 
of the machine voltage in closed loop can be defined as 
follows : 

5 Equation 3 

= a[v t (t) ,S] + b- v„W 

where a[v t ,5] is a function derived after appropriate 
substitutions and simplifications. 

The variable v RJF is the voltage reference signal for 

10 the regulator AVR such that an adequate processing of the 
error e = _ v t will allow to zero it over time, 

generally asymptotically. This form of the machine voltage 
equation is appropriate for the synthesis of the 
stabilising compensator PSS according to the invention. 

15 As a first instance, we will suppose that the 

measurement of the derivative of the load angle (d<5/dt) is 
available; i.e. the difference between rotor velocity and 
grid frequency. This allows to simplify and to understand, 
in relatively simple fashion, the functionality of the 

20 innovative technique. Thereafter with a trans format ion of 
variables, it will be indicated how it is possible to 
obtain, practically, the same results with measurements 
that are definitely more readily available on real plants. 
To complete the model that allows to perform the 

25 synthesis in subject, it is necessary to define the 
balance equation of the shaft powers and the relationship 
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20 



between electrical power and load angle. 

Given that: 
T M = ±— SL e a>(t) = $>-.S + a> R 

in which J is the moment of inertia of the axis line, a> 0 
is nominal pulsation, P 0 is nominal power and 8 is the 
nominal angle, the balance equation is: 

Equation 4 



. ,5 + a) 



the mechanical power supplied 
P B (t) = — • sin ((5 0 • 8) the electrical power generated 



X B 



Replacing the latter expression in the balance 
equation, an alternative relationship for machine voltage 
is obtained: 



Equation 15 

v t (t) = v t (S) = 



*>0 



sin(<y o • 8) 



This allows to express machine voltage as a function 
15 of the angle, provided that 8*0; this artifice will 
allow, as we shall see, to describe the system under study 
in the canon form, necessary for the synthesis of the 
control . 

The technique according to the invention comprises 
the definition of a, controller whose objective is to zero 
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a variety (surface or multi-plane) according to a 
principle called "Sliding-Modes" . This surface is nothing 
other than a linear combination of the states that 
describe the dynamics of the system, provided that said 
5 states are represented in the form: 

Equation 6 

Xj = x 2 

X 2 ^3 

x n = /(x) + g(x) • u 

X6i? w 

the canon form. 

If f (x) is uncertain but limited and g(x) is 
10 different from zero and has known sign, a control: 

Equation 7 

u{t) = -A:(x) 2 • sign(cr) and 

<?(t) = I-t- + A x x = c x • x, + c 2 • x 2 + ... + x. 



( d V 



n 



and its derivative 
d«) = grad(<r) • i, 

if fcfc/ 2 is as large as to comply with the relationship 
\ d<y 2 , , 

2 " ~dT < ~ 7 * W 

15 is such that we have <r -> 0 and cr — > 0 in a finite time. 

In these conditions there is an attractive motion on 

a = 0 which is called "sliding-modes" . 

Therefore, the system is defined in the required 

form: 
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Equation 8 

dS A 

~dT = 5 

d A = 8 
dt 

~ = F[8) + G{S) -v^ 
in which: 



Equation 9 

P M - — • [fl(v, f J) • sin(£) + <? 0 ■ <j • v, • cos{<y)] - T u • — • S* 

o 0 



F(S) = 



• sin(<5) 

G(<?) = - — ^ . fe(v, , S) 



£■'♦•0 



having already defined v t in equation 5 and v t in equation 
3. 

Equation 8 represents a system which, in fact 
(complete model) is, under normal operating conditions, 
locally stable though with oscillating phenomena; this 
assures that if the control generates an attractor on 
o=0, then motion according to "sliding-modes" will be 
obtained even if k 2 meets the given specifications only 
around o=0. 

The surface selected for the compensator PSS is the 
following: 
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Equation 10 

cr(t) 



= ^ + • 8 , X > 0, which yields 



&{t) = + A . = W) + W) . v + 2 . ^ 



The substitution of the control law in Vrxf determines the 
following equation: 

5 Equation 11 

<x = F(S) - G{S) • * 2 • 5igw(<r) 

which, if k is chosen appropriately, yields: 

1 da 2 . , 
a <> -7} - sign{cr) => — • < -77 • \<y\ 

2 dt 

which is the condition required by S.M. 

10 In the present invention, the same principle is 

implemented, but making use of a particular description of 
state variables and, moreover, the control law is in a 
form that allows, at the designer's discretion, to obtain 
or not a dampening of the oscillations which are typical 

15 of this technique. 

Since the measurements of network frequency are 
unlikely to be available, the following position: 

Equation 12 

P V 
% = iJL = 12L . S i n ( S) 

v t x E 

20 allows to reach a similar form to the one obtained in Eqn. 
8; in fact 
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Equation 13 

y 

£ = — • cos ($) S 0 • 8 

which, in the normal operating range, is: 
4 = S 0 - - ^ • ij ■ da cui si ha S. 

Operating successive derivations and substitutions, 
the following equation system is reached: 

Equation 14 

dt h 
*i = p 

dt * 

^ = <J><|, 1 4) + T(4. 1 4) • 

equivalently to the system 8 . 

If one selects in this case, too, a surface: 

Equation 15 

<r{t) = + JL^j • £(t) we have, as before : 

<Ht) = ^f + A- 4 = ® + r- VRlF +A-4 
dt 

Obtaining the elimination of this surface and of its 
derivative, entails: 

Equation 16 

dt * 

-d7 = <y - x -t 

dt 
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which, by definition of means also eliminating d&dt and 
d 2 #dt 2 as desired; in fact: 

Equation 17 

V 

€ = 0 => # = constant = — - sin (8) => 8 = constant 
in quanto, dalle ipotesi fatte, e x F sono costanti. 

in so far as, form the hypothesis made, V R and X E are 
constant . 

Therefore, we have: 

Equation 18 

8 = 0 and 
8 = 0 

asymptotically, therefore, from equation 4, we have: 

The realisation of the present invention requires a 
synthesis of a de facto estimated surface since the 
available measurements are not sufficient for its 
construction. 

The techniques whereon the design of the invention is 
based are essentially aimed at estimating the sliding 
surface. 

First of all, we observe that during the normal 
operation of an alternator, the angle 8 is expected to be 
limited in a range between about +/- x/6 and the angular 
frequency of the machine always very close to the nominal 
value, it is possible further to simplify the initial 
equations : 
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Equation 19 



d 2 S ^ P M - P B 



u2 

V 



dt 2 T M 



P E » — j- • v t • 8 where 



- S 0 ^ 5 <* S 0 con S 0 & and xl = ^ 



6 ~ "* sin(£ 0 ) 

These positions allow to manipulate equations that 
are definitely simpler, without thereby losing in 
generality. 
Therefore: 

Equation 20 

Vr Vr 



20.1) — |- = j-Z da cui v ( . Derivando si ha 

dt 2 x 1 



20.2 ) d -4- = - ' Vf + — - avendo posto /» « 0 e, infine 
dt 3 T. A 



v, 



20.3 ) 



Now remembering the equation of the machine voltage 
10 (Equation 3) and replacing in 20.3: 



Equation 21 



V R 



in which b>0. 
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As is readily apparent, since all quantities aire 
uncertain (V R/ x' E /T M / etc) in the so-called "drift" term, 
(i.e. the first term at second member of the equation), 
the trans format ion of variables made allows to apply the 
5 "sliding modes' 7 technique which is practically insensitive 
to the variation and/or uncertain of these parameters. 

The synthesis of the sliding surface requires to 
estimate two system states: the first and second 
derivative of the variable 
10 Figure 2 schematically shows the plant 100 

highlighting the implemented regulation architecture. 

Some of the blocks shown in Figure 2 represent 
components shown in Figure 1 and, therefore, they were 
designated with the same symbolic reference. The block 
15 designated with the symbol GEN & NET schematically 
represents the alternator G and the network NET defined 
above. The block designated as PAR represents the 
parameters of the network NET (network voltage V R/ 
external reactance Xe, and grid frequency F R ) . Figure 2 
20 also highlights the flow of the signals between the 
various blocks and it is in accordance with what is 
described with reference to Figure 1 and with what is 
expressed in the analysis set out above. 

Figure 3 shows a particular embodiment of a 
25 stabiliser device PSS according to the invention and 
usable in the plant 100 of Figure 1. Preferably, the 
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stabiliser device PSS is embodied as a modular element by 
means of analogue and/or digital components housed on a 
board 15; in particular, the stabiliser device PSS can be 
selectively coupled to and removable from the voltage 
regulator AVR and is thus interchangeable to be replaced, 
for instance in case of failure or for upgrading 
operations . 

According to this first embodiment, a second order 
deriving filter is used. 

The stabiliser device PSS is provided with a first 1 
and a second 2 input terminal for its main input signal 
which, according to said example, are active electrical 
power P E and machine voltage v t . 

The stabiliser device PSS implements the following 
equations: 

Equations 22 



Out _ PSS = H • signg) • /(&) • sign(a) 

The equations 22 are similar to the equations 15 and 
to the first of the equations 7. In particular, the symbol 
x>/v " indicates that the corresponding quantities are 
obtained by estimation. Moreover, it should be recalled 
that the estimated quantity £ is expressed in accordance 
with the equation 12 . 

In the second of the equations 22, note that the gain 
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H has the sign of £ to guarantee the sliding modes 
condition and f (a) is any suitable function. 

For example, for the sliding modes approximation the 
function f(a) is: 
5 f=l, |a|/flfe>l;f - Icr|, |ct|/0<1. 

A function of this kind is indicated by the 
aforementioned text by Slotine, "Applied Nonlinear 
Control" and it is such as to allow advantageously to 
dampen "chattering" . 
10 Returning to Figure 3, the input terminals 1 and 2 

are connected to a divider block 3 which is such as to 
divide the signal P E and v t and to obtain the estimate of 
the variable £(t) in accordance with the equation 12. 

One output of the divider block 3 is connected to 
15 first processing means 4 which conduct a processing 
operation to provide an electrical signal indicating an 
estimate of the sliding surface cr (t) . 

According to the first embodiment of the invention, 
the first processing means 4 include a second order filter 
20 able to conduct a processing operation in accordance with 
the first expression of equation 22 which includes an 
operation of computing the first derivative of £(t) and 
the second derivative of the same variable £(t) . 

Figure 3 shows, purely by way of example, a possible 
25 transfer function of the second order filter 4. In 
particular, the filter 4 can be implemented by an analogue 
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or digital high pass filter . 

An output OUT of the second order filter 4 is 
connected to a processing block 5 such as to evaluate the 
sign of cr(t) . 

5 Moreover, the output OUT is connected to a processing 

branch including a plurality of blocks 6-10. 

The block 6 is such as to evaluate the absolute value 

of the signal indicative of the sliding surface cr(t) and 

is connected to an additional divider block 8 such as to 
10 effect the ratio between the output signal of the block 6 

and the parameter 0 introduced above (block 7) , set 

according to the example to the value 0.01. 

The divider block 8 is connected to a saturation 

block 9 which saturates to 1 the ratio <y(t)/3>. 
15 The output of the saturation block 9 is connected to 

a gain block 10 which introduces the aforementioned gain 

H, for instance, set to the value 0.5. 

The output of the gain block 10 and the output of the 

block 5 are connected to a first multiplier 11 having an 
20 output which in turn is connected to a second multiplier 

12. 

The second multiplier 12 is such as to receive at its 
input also a signal corresponding to the sign of the 
estimate of the quantity £(t) made available by a block 13 
25 connected to the output of the divider block 3. 

Moreover, said second multiplier 12 allows to 
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multiply the signal sign(£(t)) with the signal H f(<r(t)) 
sign(cr(t)) output by the first multiplier block 11 and to 
generate the output signal OUT_PSS present on the output 
terminal 14 . 

5 In the operation of the compensator device PSS, based 

on the signal indicative of the active power P E and of the 
one for machine voltage v t/ the sliding surface cr(t) , of 
the second order sliding modes type, is estimated by the 
block 3 and by the filter 4. Moreover, the block 13 

10 evaluates the sign of the ration P E /v t , i.e. the sign of 
the quantity £(t) that constitutes the variable to be 
controlled. 

Starting from the signal corresponding to this 
estimated sliding surface a(t) the function f(cr(t)) is 

15 built (block 6-9) and the sign of the signal corresponding 
to said surface is evaluated. 

Through a first multiplication in the first 
multiplier block 11 and a subsequent multiplication in the 
multiplier block 12, the signal OUT_PSS is built. 

20 This signal OUTJPSS, supplied to the summing node ND 

of the exciter 200 of Figure 1, will contribute to the 
formation of the voltage reference signal V RIF which is 
inserted in the voltage regulator AVR. The action of the 
voltage regulator ACR and of the actuating device ACT 

25 operatively associated with the field windings CL-F causes 
an intervention on the operation of the alternator G. 
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In this way the measuring signals P E and v t will be 
modified causing the amplitude of the signal £(t) obtained 
in the stabiliser PSS to be maintained substantially 
constant (i.e. the first derivative of £(t) is 
5 substantially equal to zero) and the signal a(t) to be 
brought substantially to converge to zero. 

This allows to stabilise the unit 50 damping its 
electromechanical oscillations. 

According to a different version of the embodiment of 
10 Figure 3, the stabiliser PSS is also supplied with a 
signal indicative of the required active power Prif and a 
signal indicative of the required machine voltage V SE t- 

In this- case, the sliding surface cr(t) will contain 
an additional term which will be directly constituted by 
15 the variable, £(t)=£(t)= P E /v t , i.e.: 

Equation 23 

a(t) = (j| + aJ • fan - fiurWl 

in which ft^tf) = Pjup/^set and P RXF and V SET are, respectively, 
the reference values of active power and machine voltage 
20 desired at steady state. 

This last assumption allows to have an easily 
computed surface: 

Equation 24 

<*> = [j. + A J • [ft) - W>] = ^ + 2 • A • *8j± + A 2 . [ft) - W>] 
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The rest of the control is similar to the one 
described in Figure 3 . 

Figure 4 shows a second embodiment of a PSS 
compensator in accordance with the present invention. 
5 According to the embodiment of Figure 4, the processing of 
the measurement signals for estimating the sliding surface 
cr(t) is performed not with a second order filter but with 
a linear or non linear observer 4' . Moreover, it should be 
noted that the observer 4' can be obtained with a sliding 
10 modes technique (sliding observers) • 

In the diagram of Figure 4, the same numerical 
references of Figure 3 are used to designate identical or 
similar components . 

The observer 4' is connected to the output of the 
15 divider block 3 to receive the signal P E /v t and to a third 
input terminal 2' to receive at least another signal. 

For example, the third terminal is provided with the 
voltage reference electrical signal v RIF which is also sent 
to the voltage regulator AVR shown in Figure 1. Note that 
20 the fact that to the compensator PSS is sent the reference 
signal of v R i F is also indicated in the diagram of Figure 
2. 

The synthesis of a non linear observer exemplifier 
can easily be made by a person skilled in the art, based 
25 on the following expression: 
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Equation 25 

£ = 6 

I, = 4 + a • « - &) 
-A = & + • m - &) 

which, once gi, g 2 e g 3 , are properly selected, will allow 
to eliminate the observation error and will provide an 
5 estimate of the unmeasured states. 

Moreover, it should be noted that the stabiliser PSS 
can be obtained, according to a variant of the invention, 
employing instead of the filter 4 or of the observer 4' a 
Levant differentiator, i.e. a device also known as "super- 
10 twisting 2-Sliding modes dif f erentiator" . 

Figure 5 shows a compensator PSS according to a third 
embodiment of the invention, in which control is based on 
the measurement of rotor velocity a>. This embodiment is 
particularly advantageous because the measurement of rotor 
15 velocity co is commonly available in the velocity 
regulators of turbo-machines and hence can be easily 
acquired in rapid and precise fashion. 

In this case, one starts from the equations in 8: 
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Equation 26 

^- = a ) -co R = Sif) 
at 

— = Cm - C " con C, = — • (y a ■ cos 5 + y/ d • sin 5) 
dt J x. 



d 2 S 
dt 2 



C m - ^ • jV, • cos 5 - • sin 5) • <* + • sin 5 + ^ • cos /J 



where, as always, J is the angle between axis q and 
5 the network voltage vector, having amplitude V R , a> is the 
angular velocity of the reference system integral with 
axis q (i.e. rotor velocity), cfo is electrical network 
frequency, and lastly C m and C e designate the motive torque 
and the electrical torque, respectively. Moreover, y/d and 
10 y/q are the direct axis flow, and the quadrature axis flow, 
respectively, and they are defined as follows: 

Equation 27 

dy/ d 1 ,1 x d di d 

dt v d0 Vi r d0 f r d0 dt 
dt 7\o ^ + r q0 '< x * dt 

co • y/ d — v R • cos 5 
v R • sin £ + fi> • y 9 

where i d is the current associated with the axis d, i q is 
15 the current associated with the axis q, x e/ x d , x d ' , x q/ 
x q ' are, respectively, the external reactance, synchronous 
direct axis reactance, transient direct axis reactance, 
synchronous quadrature reactance and transient quadrature 



WO 2005/013462 



PCT/IB2004/002480 



30 



reactance, T' d0 , T' q0 the load-less direct axis and 
quadrature time constants and v f is field voltage, which 
represents the control variable. 

The equations 26 can be manipulated so as to express 
5 the field voltage v f and to identify the possibility of 
obtaining the convergence to zero of the selected sliding 
surface. For this purpose, some simplifying hypothesis are 
made, which leave the generality of the approach unchanged 
and are commonly verified. 
10 In the first place, under normal operating conditions 

rotor frequency co is about equal to network frequency (Dr. 
Moreover, the disturbance represented by mechanical torque 
can be considered constant. The first and the third of 
the equations 26 therefore assume the following form: 

15 Equation 28 

8 = co - co R « 0 



The quadrature flux y/ q relating to the axis q is 
represented by a differential equation of relative order 
zero which can be replaced, with an acceptable 



Therefore, the time derivative of the quadrature flux yr q 
is proportional dS/dt which, as shown by the first of the 
equations 28, is about nil. 

Consequently, the following apply: 




dt 2 



J 



20 approximation, by an algebraic equation in sin 5. 
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dt 7 



Equation 29 



4 d>Va 

• sin 8 + — r 2 - • cos 8 



dt 



dS 



dO 



V, + 



7 *-r • V R • COS 8 



dy/ q 



- sin 5 and 



■ sin £ 



dt ~ [x e + x q • a) ' " 

At this point we can derive an expression of the 
direct axis flux y d recalling that: 



Equation 30 



d8_ = C m - 



C 8 _ 



C m ~ — • (y g • cos <? + • sin <?) 



Vd - 



cos £ 



sin £ 



wherefrom the following is derived 



<e • ( X e + X ' &) 



C„ - J 



V* = 



dt 



+ V R ' X q ' S ^ n & ' COS & 



v R • {x e + x q ■ co) • sin S 

Therefore, the direct axis flux y/ d , the quadrature 
flux y/ q and the derivative of the quadrature flux yr q are 
functions of the angle. 8, whereas the derivative of the 
direct axis flux yr d is a function both of the angle 5, and 
of the field voltage v f/ i.e.: 
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Equation 31 

Wa = WAS) 

= a d (S) + b d ■ v f 



d W* _ 



dt 

Replacing in the Equation 29, the following is 
obtained: 

Equation 32 

_ . ( a<f (£) . S in 5 + a tf <£) • cos S + 6 rf • sin 5 • v r ) 

fL* « _3 = F(<5) + g(<J) • v y 

<ft 2 «/ 



with 



F(<5) = - ^ 



Vp a,(£) • sin £ 4- a^(^) ■ cos 8 



g(S) = 



b d » sin J 



In normal operating conditions, the drift term F f# 
10 is a limited quantity and g (S) has negative sign for 
0<8<n/6 (i.e. still in normal operating conditions). 
Consequently, the conditions for the sliding motions 
exist . 

A sliding surface and the relative control law 
15 (expression of the field voltage) are given by the 
following: 

Equation 33 

, dS , . 

a* = A * « A • €0 

dt 

v f = A 2 - s/gw(tf>) (since g (5) < 0) 

with the usual meanings of the symbols. 
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The type of control described, based on rotor 
velocity <o (or, more precisely, on its derivative a>) is 
accomplished by the compensator PSS of Figure 5. In 
detail, in the third embodiment, the compensator PSS has 
5 an input terminal 20 connected to the measurement 
acquisition and processing module ACQ-M to receive the 
rotor velocity a). A differentiator block 21 receive the 
rotor velocity a> through the first terminal 20, determines 
the derivative a> of the rotor velocity a) and supplies it 
10 to a gain block 22, arranged immediately downstream, which 
operates a multiplication times the factor A. In practice, 
the differentiator block 21 and the gain block 22 together 
calculate the sliding surface a based on equation 33 . The 
sliding surface a is supplied to a processing branch 24 
15 which includes a divider block 25, an absolute value 
extractor block 26, a saturation block 27 and a multiplier 
block 28, mutually connected in cascade in the above 
order. In particular, a first input and a second input of 
the divider block 25 receive the sliding surface a and, 
20 respectively, a parameter 0 which lies within a 
predetermined range. On the output of the divider block is 
therefore present the value o/<P, the module \cr/0\ 
whereof, possibly limited by the saturation block 27, is 
calculated by the absolute value extractor block 26. The 
25 multiplier block 28, lastly, multiples the value provided 
by the saturation block. 27 for a positive gain H 2 . An 
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additional multiplier block 30 has a first input connected 
to the output of the multiplier block 28 and a second 
input connected to the output of the deriving block 21 , 
through a sign recogniser block 31, which determines the 
sign of the sliding surface. The output of the multiplier 
block 30 forms the output terminal of the compensator PSS 
and provides the output signal OUT_PSS in accordance with 
the equations 33. 

Alternatively, the derivative of the rotor velocity a> 
can be determined using a second order sliding modes 
observer (" super- twisting 2-Sliding modes differentiator 7 ' 
o Levant differentiator) . 

According to a fourth embodiment of the invention, 
the PSS compensator uses an estimate of the derivative of 
the angle 8, which is provided by the following: 

Equation 34 
8 = a> - A 3 

where A is a constant parameter. In this case, the 
corresponding sliding surface a is: 
Equation 35 



If rotor velocity a> us constant, the estimate 5 
tends to zero at a rate that depends on the parameter. It 
has already been shown that, if the control law is such as 
to reduce to zero the sliding surface a and its 
derivative, the following applies: 
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Equation 36 

d 2 8 . 3 d8 
dt 2 dt 

Moreover, the system is autonomous and asymptotically 
stable. When, instead of the derivative 8, its estimate 8 
is used, we have: 

Equation 37 

d 2 8 . ^ dS 



= & _ U - A) . ™L + (a - i 



* 2 ""' ^ 

from which it may be deduced that it necessary to choose 
0<K<\. 

10 The control described above is implemented by the 

compensator PSS of Figure 6, which comprises the 
processing branch 24, the multiplier block 30 and the sign 
recogniser block 31, already described with reference to 
Figure 5, and further includes an input terminal 40, a 

15 subtr actor node 41, an integrator block 42, a gain block 
43 and a processing block 44. In detail, the subtractor 
node 41 has a first input connected to the input terminal 
40, which receives the measurement of the rotor velocity a> 
from the measurement acquisition and processing module 

20 ACQ-M, and a second input connected to the output of the 
gain block 43. The integrator block 42 has an input and an 
output, respectively connected to the output of the 
subtractor node 41 and to the input of the gain block 43, 
so as to form a loop. In practice, the estimate 8 of the 

25 derivative 8 of the angle 8 is provided on the output of 
the subtractor node 41, and is used by the processing 
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block 44 to compute the sliding surface a in accordance 
with equation 35. The sliding surface a is then processed 
by the processing branch 24 and by the multiplier block 30 
as described above with reference to Figure 5. 

A fifth embodiment of the invention provides for the 
use of the active electrical power P E/ which, like rotor 
velocity a>, is easy to measure and is normally available 
in electrical power production plants. The approach 
described below can be advantageously exploited under some 
simplifying hypotheses, which are generally verified. In 
particular, it is assumed that the network parameters, 
with the exception of the load angle 8, are constant or 
slowly variable and that the active power P E depends on 
machine voltage v t , on network voltage V R and on the angle 
5. Consequently, the active power P E is also a function of 
the direct axis flux y/ df which in turn has a first order 
dependence on the field voltage v f . In fact, the active 
power P E is also correlated to the quadrature flux \j/ q , 
which can, however, be expressed in terms of the angle S. 
Therefore, we have, in practice: 

Equation 38 

P E = ¥ d ) 

Lastly, it is assumed that the following 
relationships apply: 
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Equation 39 



S » 0 



5 m ?M = Pg 

The sliding surface a can be of the general type 
described by the following: 

Equation 40 

or = 1 - P F 

which, taking into account the equations 38 and 39, 
entails : 



da 


X • — 


~dP s 


8 + 


dP B 


■ V^j = 




dt 


dt 


38 






= X • 


~d 2 P E 
dS 2 


8 2 + 


dP E 

dS 


■ 8 + 


d 2 P E . 2 

- 2 * V<* 





and: 



Equation 41 

r 



da 
dt 



= 2 



d 2 P B 



dP„ 



dP E 



The control has a 1 st order dependence on the direct 
axis flux y d ; therefore, the control law, which is 
represented herein again by the field voltage v f , is of 
the type: 

Equation 42 

v f = J - k 2 • sign(a) • dz 



The control term is therefore the following: 
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Equation 43 

dv f 

—f- = -k~- sign[a) 
dt 

The control described herein is implemented by the 

compensator PSS of Figure 7, the structure whereof is 

5 similar to the device shown in Figure 5. In particular, 

the compensator PSS in accordance with the fifth form of 

embodiment of the invention comprises the differentiator 

block 21, the gain block 22, the processing branch 24, the 

multiplier block 30 and the sign recogniser block 31. An 

10 input terminal 50 is connected to the measurement 

acquisition and processing module ACQ-M, to receive the 

measurement of the active power P E . Moreover, an 

integrator 51 and a saturation block 52, which provides 

the output signal OUT_PSS, are connected downstream of the 

15 multiplier block 30. In practice, the differentiator block 

21 and the gain block 22 determine the sliding surface <x 

on the basis of the equation 40 and feed it to the 

processing branch 24 and to the sign recogniser block 31. 

In accordance with the equation 43, therefore, on the 

20 output of the multiplier block 30 there is present the 

value — — = -k 2 • sign{cr) , wherein the sliding surface a is 
dt 

the one calculated by the deriving block 21 and by the 
gain block 22. Lastly, the integrator block 51 calculates 
the field voltage v f to be provided at the output and the 
25 saturation block 52 limits the extent of the control, if a 
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pre-determined threshold is exceeded. 

According to a sixth embodiment of the invention, 

whereto Figure 8 refers, the compensator PSS uses an 
estimate of the derivative of the active power P E . For 
this purpose, a second order filter is used which allows 
to construct a suitable sliding surface cr, with no need to 
perform derivatives of the measurements and to integrate 
the control signal. In particular, the sliding surface cr 
is given by the following: 

Equation 44 



where the estimate —— of the derivative of the 



active power P E is calculated by means of a second order 
filter 60 (see Figure 8) which uses the following: 
Equation 45 



In particular, the parameters a and b of the filter 
are chosen in such a way that the estimate error 





e = P E - P E asymptotically tends to zero, i.e.: 



Equation 46 
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If the parameters a and b of the filter allow to meet 
this condition, it is possible to construct a 
discontinuous and canonical first order control directly 
on the field voltage v f (instead of on its derivative) . 

5 As shown in Figure 8, in this embodiment of the 

invention, the filter 60 replaces the differentiator block 
21 and the gain block 22 of Figure 7* Moreover, the 
integrator block 51 and the saturation block 52 are 
missing: in this case, the control variable, i.e. the 

10 field voltage v f , is provided directly by the multiplier 
block 30. 

Figure 9 shows a seventh embodiment of the invention, 
whereby the compensator PSS uses, in addition to the 
active power P E/ also an estimate or a measurement of 
15 mechanical power P M or, alternatively, an estimate or a 
measurement of the mechanical torque C M . The compensator 
PSS of Figure 9 comprises, as in the previous cases, the 
processing branch 24, the multiplier block 30 and the sign 
recogniser block 31. Moreover, upstream of the processing 
20 branch 24 is positioned a module 7 0 for calculating the 
sliding surface o> which has a first and a second input 
terminal 71, 72 connected to the measurement acquisition 
and processing module ACQ-M, to receive the rotor velocity 
a) and the active power P E . In detail, the module 70 
25 comprises an observer 73, a calculation block 74, a 
multiplier 75 and a selector 76. 
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The observer 73 provides an estimate of the balance 
between the electrical torque C E and the mechanical torque 
C M/ based on the rotor velocity co and on the active power 
P E . The calculation block 74 has an additional input 77 
5 connected to the measurement acquisition and processing 
module ACQ-M, to receive the measurement of the mechanical 
torque C M and, still based on the rotor velocity CD and on 
the active power P B , it calculates the balance between the 
measurements of the electrical torque C E and of the 
10 mechanical torque C M . 

The selector 7 6 alternatively connects the outputs of 
the observer 73 and of the calculation block 74 to the 
multiplier block 75, which multiplies the quantity 
received at its input times the parameter X. 
15 Preferably, the observer 76 is in turn provided with 

a sliding modes control and it estimates the mechanical 
power P M as if it were a constant (the hypothesis is 
realistic, since the mechanical power P M is slowly 
variable) „ 

20 Starting from the equation of motion 

Equation 47 

da_ = C M - C E 
dt J 

the observer 73 is constructed as follows: 
Equation 48 

25 = ^ = C * + G -e 

dt J' 

in which the moment of inertia J' is unknown, the 
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electrical torque C B is measured, G is the observer gain 
and e = co — & is the error of the observer 73. 

The sliding surface obss of the observer is given by 
the following: 

Equation 49 



■f- 



^ 

e 



whereas the control law of the observer is given by the 
Equation 50 

dt 

The error of the observer tends to zero, as will be 
shown, and therefore the difference between the real 
torque balance and the estimated torque balance tends to 
be reduced to zero, thus asymptotically one has: 

Equation 51 

J' J 

Starting from the estimate error, we have: 
Equation 52 

e = a> - & = — — — — — — g • e 

J %J* 

from which: 



Equation 53 

C C 



fx x\ 
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and, replacing the equation 50 in the equation 53: 



Equation 54 

8 = - ^1 - . C E - g • e - ^ . rigger) = d 05S - ^ 



*oss 

Hence, since the conditions for the existence of a 
.5 majorant as required for the sliding modes control are 
met, we will have the motion on the sliding surface cross, 
which tends to zero together with its derivative & oss . 

Therefore, based on the equation 52, the error e of the 
observer 73 also tends to zero asymptotically and the 
10 equation 51 is obtained, which assures the correctness of 
the estimate made by the observer 73. 

The compensator device PSS of the invention has 
considerable advantages. 

One of these advantages is that the compensator 
15 device PSS of the invention does not require, for its 
implementation, the full knowledge of the process 
parameters to be controlled. 

As can be inferred by the above description, the 
output signal OUT_PSS of the compensator PSS is obtained 
20 based on functions linked to the sign of the variable to 
be controlled £(t) and to its derivatives and not based on 
functions which are directly linked to the values assumed 
by process parameters. 

Unlike traditional compensators, this assures a high 
25 robustness and tolerance of the compensators according to 
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the teachings of the invention with respect to 
uncertainties on the process parameters. 

Moreover, the choice of a first order sliding modes 
control is particularly advantageous in terms of 
5 simplicity of construction and robustness of the 
compensator . 

It should also be observed that the examples of 
embodiments of the invention described above are based on 
two additional signals to be fed to the compensator PSS 

10 (e.g., active power P E and machine voltage v t ) already 
available on every energy production plant, making the 
implementation of the teachings of the invention 
particularly easy* 

Note also that it has been noted that the undesired 

15 oscillations of the output signal of the compensator 
0UT_PSS ("chattering" phenomenon) , intrinsically present 
in the sliding modes control circuits, are not relevant in 
the context of energy production plants. It has been 
observed that these oscillations (at relatively high 

20 frequency) are filtered by the actuators of the exciter 
200 and hence do not negatively influence the desired 
regulation. 



